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Abstract. We present an overview of the upcoming Microwave Anisotropy Probe 
(MAP) mission, with an emphasis on those aspects of the mission that simplify the 
data analysis. The method used to make sky maps from the differential tempera- 
ture data is reviewed and we present some of the noise properties expected from 
these maps. An overview of the method we plan to use to mine the angular power 
spectrum from the mega-pixel sky maps closes the paper. 

1 Mission Overview 

In 1992 NASA's Cosmic Background Explorer (COBE) satellite made a full 
sky map of the cosmic microwave background (CMB) temperature with 7° 
resolution, uncorrelated pixel noise, minimal systematic errors, and accurate 
calibration, from which CMB temperature anisotropy was first discovered ||l|, 
07 [§1 ■ The purpose of the MAP mission is to re-map the anisotropy over 
the full sky with more than 30 times the angular resolution (~ 0.23° FWHM) 
and more than 35 times the sensitivity ('-^ 20 fiK per 0.3° pixel) of COBE, 
but with the same level of quality control as was possible with COBE. With 
this data, MAP will measure the physical interactions of the photon-baryon 
fluid (sound waves) in the early universe and thereby test models of structure 
formation, the geometry of the universe, and inflation. 

In the years since COBE, a host of ground-based and balloon-borne exper- 
iments have detected and characterized fluctuations at smaller angular scales, 
most recently the experiments TOCO §, BOOMERanG || and MAXIMA 
[0. However, because of their proximity to the Earth and its atmosphere, 
none of the ground or balloon-based experiments enjoy the extent of system- 
atic error rejection or calibration accuracy that was possible with COBE. 
Moreover, many of these experiments have significantly correlated noise that 
places severe demands on the data analysis. The COBE data still serves as 
a benchmark for the field and many aspects of the COBE mission have in- 
fluenced the design of the MAP mission. The need to minimize the level of 
systematic errors has been the major driver of the MAP design and has led 
to the following high level design features: 

• A highly symmetric diflferential design: MAP is a differential exper- 
iment based on pseudo-correlation microwave radiometers that employ 
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phase-matched HEMT amphfiers. The instrument measures temperature 
differences between two points ^ 141° apart on the sky. By measuring 
temperature differences, rather than absolute temperatures, many spuri- 
ous signals will be common-mode and thus cancel upon differencing. Also, 
by employing a pseudo-correlation design with a fast chopping frequency 
between the two sky inputs, 1/f noise that arises from the HEMTs can 
be chopped out. The resulting power spectrum of the radiometer noise is 
very nearly white (see §3). 

• Multi-frequency: There are five frequency bands from 22-90 GHz that 
will allow emission from the Galaxy and other non-cosmological sources 
to be modeled and removed based on their frequency dependence. In the 
lowest frequency bands, MAP will probe the high frequency tail of radio 
emission from our Galaxy and provide valuable data on the enigmatic 
microwave foreground emission that correlates with thermal dust emis- 
sion, but has a much different spectrum. See ^ for a recent summary of 
evidence for this foreground. 

• Stability: MAP will observe from a Lissajous orbit about the L2 La- 
grange point 1.5 million km from Earth. The L2 point offers an excep- 
tionally stable environment and an unobstructed view of deep space, with 
the Sun, Earth, and Moon always in shadow behind MAP's Sun shield. 
map's large distance from Earth protects it from near-Earth emission 
and other disturbances. While observing at L2, MAP's Sun shield and 
solar panels maintain a fixed angle with respect to the Sun to provide 
exceptional thermal and power stability. 

• Low beam sidelobe levels: The MAP optical system was designed with 
the foremost goal of providing adequate angular resolution along the line 
of sight while at the same time rejecting stray light from other directions. 
For example, the largest instantaneous signal due to radiation from the 
galactic plane spilling into a sidelobe is expected to be less than 2 /iK at 
90 GHz. 

An overview of the MAP satellite is shown in Fig. |l]. The major visible 
features of MAP include the back-to-back telescope optics with 1.4 x 1.6 m 
primary mirrors and 1 m secondary mirrors, the passive thermal radiators 
which cool the HEMT amplifiers to <100 K, the hexagonal structure housing 
the spacecraft service modules, and the large solar panel array/Sun shield 
which keeps the instrument in full shade. MAP weighs a total of 830 kg 
and stands ~4 m tall. It will be launched in 2001 aboard a Delta 7425-10 
expendable launch vehicle from the NASA Kennedy Space Center Eastern 
Test Range. 

2 Scan Strategy 

The MAP scan strategy plays an important role in systematic error rejection. 
It was designed with the following goals in mind: 
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Fig. 1. An overview of the MAP satellite. 



• Scan a large fraction of the sky as rapidly as possible, consistent with 
reasonable requirements on the controlling hardware and the telemetry 
data rate. 

• Scan each sky pixel through as many azimuthal angles as possible for the 
reasons listed below. 

• Observe a given pixel on as many different time scales as possible. 

• Maintain the instrument in continuous shadow for optimal passive cooling 
and avoidance of stray signals from the Sun, Earth, and Moon. 

• Maintain a constant angle between the Sun and the plane of the solar 
panels for thermal and power stability. 

The strategy that was ultimately adopted combines a "fast" spin about the 
spacecraft symmetry axis with a slow precession about the Sun-MAP line 
(which is always within 0.1° of the Sun-Earth line at L2). Since each telescope 
line of sight is ~ 70° off the symmetry axis, the path swept out on the sky 
by a given line of sight resembles a spirograph pattern that reaches from the 
north to south ecliptic poles. As MAP orbits the Sun, this pattern revolves 
around the sky so that full sky coverage is first achieved after 6 months of 
observing at L2. 

The MAP scan strategy achieves a reasonable level of azimuthal coverage 
in each sky pixel. For example, a pixel in the ecliptic equator is observed 
over ~30% of the possible angles of attack; a pixel at the cusp of the annular 
coverage at ~ 45° ecliptic latitude is observed over about 70% of possible 
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Fig. 2. A schematic of the MAP scan strategy from L2. The satellite covers approx- 
imately 35% of the sky each day. The scan pattern on the right depicts the motion 
of a single line of sight over a 1 hour precession period (lighter spirograph pattern), 
and a ~2 minute spin period (darker single circle). The inner circle depicts the path 
of the spin axis during one precession. 



angles of attack; and a pixel near the ecliptic poles is observed from 100% 
of the possible azimuthal orientations. A large azimuthal coverage provides 
numerous desirable features in the data: 

• Helps to produce a stable sky map solution. 

• Produces small pixel-pixel covariance at the beam separation scale. 

• Minimizes striping due to any residual 1/f noise in the differential data. 

• Maximizes polarization sensitivity. 

• Maximizes azimuthal symmetry of the beam response on the sky. 

Prom this standpoint, the MAP strategy is not as complete as COBE's, which 
achieved nearly 100% azimuthal coverage in all pixels. However, in order to 
achieve such completeness, the spacecraft spin axis must ultimately point to 
every pixel on the sky which is less desirable from the standpoint of systematic 
error avoidance. The MAP strategy achieves reasonable azimuthal coverage 
consistent with strong systematic error constraints. Extensive simulations of 
the map making and power spectrum estimation procedures have shown that 
the strategy is more than adequate to meet MAP's scientific goals. 

3 Map Making with Differential Data 

Sky maps, t, are obtained from the differential data, d, by linear least squares 
fitting. The raw data takes the form d = A • t where A is the scan matrix 
of the experiment which has Npix columns and Nobs rows, where Npix is 
the number of sky map pixels and Nobs is the total number of differential 
observations. 
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Each row (observation) of length N^ix contains a +1 in the element (pixel) 
observed by the + horn during that observation, and a —1 in the element 
observed by the — horn. The normal equations for the sky map solution are 

(A"^ . T-i . A) • t = • • d (1) 

where T is the covariance of the instrument noise in the time domain. For 
map's differential radiometers this is well approximated by stationary, white 
noise, Tu' — {ntUt') ~ criSu'- A sample of the MAP noise covariance is shown 
in Fig. H from test data taken with the flight instrument. The whiteness of 
the noise greatly simplifies the subsequent data processing and analysis for 
two reasons 1) the normal equations for the sky map are much easier to 
solve than they would be with significant 1/f noise, 2) the pixel-pixel noise 
covariance in the final sky maps is very small and completely negligible for 
most applications. 




Fig. 3. The auto-correlation function (denoted T in the text) of 8 hours of noise 
from one of the MAP 90 GHz radiometers, taken during a cold instrument test in 
March 2000. The instrument was operating at expected flight temperature while 
observing temperature controlled targets over the flight feeds. The only processing 
applied to the raw data was to subtract a baseline with power on periods of greater 
than 1 hour (a Legendre polynomial of order 8 fit over the 8 hour period of data). 
This approximates the procedure that will be applied to the flight data. Note the 
point at lag 1 (inset) has a covariance of 1.2% of the raw noise due to a low-pass filter 
apphed just prior to the a/d conversion. There is no other significant covariance in 
the noise. 
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With stationary white noise, the maximum hkehhood sky map solution 
reduces to 

t = (A^ • A)-i • (A^ -d) (2) 
while the pixel-pixel noise covariance is 



(3) 



Formally, the sky map solution (g) 
M = (A"^ • A) which has the form 
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where the ith diagonal clement is the number of times pixel i was observed by 
either side of the instrument, and the ijth off-diagonal element is the number 
of times pixels i and j were observed simultaneously. Because of the fixed 
separation between the two instrument beams, only pairs of pixels separated 
by a fixed distance can be simultaneously observed, thus this matrix is very 
sparse. Furthermore, because the scan strategy generates ample azimuthal 
coverage, we have Nij ^ A'i so M is diagonally dominant. [Formally M is 
singular because the sum of each of its rows or columns is zero, due to the fact 
that the time-series data is differential. Thus the sky map solution and its 
pixel-pixel covariance are undetermined up to a constant since M • C = for 
any constant vector or matrix C. In practice this mode is readily projected 
out of the data in any actual inversion scheme.] 

The form of M suggests an iterative approach to solving for the map that 
employs a diagonal pre-conditioner as an approximate inverse Define 

and assume an approximate initial map solution to = t -I- 5t such that 

M • to = M • (t + 5t) = V + 5v (6) 
where v = A^ -d. This gives an expression for St in terms of known quantities 



M • (5t = M • to - V 



(7) 



We can use the above pre-conditioner to obtain an iterative estimate of the 
map correction Sti, and hence an improved solution ti 
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As was noted in jl^, an efficient way to evaluate the right hand side of 
is to group the operations as follows 



t„+i = (M^i • A^) • (d - A • t„) + 1„ 



(9) 



The interpretation of is that for each pixel the new sky map temperature 
is equal the average of all differential observations of that pixel, accounting 
for the sign of the observing horn, corrected by an estimate of the signal in 
the pair horn, based on the previous sky map iteration. The expression in 
can be efficiently evaluated because the sums can be accumulated by reading 
through the time-series data, read from disk, and accumulating data into 
arrays of length Npi^. It is never necessary to store or invert and Npi^ x Npi^ 
matrix. 
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Fig. 4. 2-point angular correlation function of a simulated MAP noise map. The 
function was computed from a residual sky map that was obtained by subtracting 
an input sky map from a solution of the map-making algorithm. The noise used in 
the simulated differential data was statistically identical to the actual radiometer 
test data shown in Fig. ^ Note the small (<1%) nearest-neighbor covariance due 
to the nearest-neighbor covariance in the time-series data. There is also a small 
(~0.1%) covariance between pixels separated by the beam spacing, indicated by 
the dashed line in the inset panel. 



The performance of this algorithm has been tested extensively with simu- 
lations that capture differential data with realistic instrument noise observed 
with the MAP scan strategy. The solution is seen to converge to the input 
in less than 50 iterations, starting from an initial guess with no anisotropy. 
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In cases where the instrument noise has been suppressed for study purposes, 
peak-peak features in the residual sky map (output — input) are less than 
0.1 /iK. In cases where realistic noise is included, the performance of the al- 
gorithm can be measured by the 2-point angular correlation function of the 
residual (noise) map. Ideally, this would be everywhere consistent with zero 
except for the bin at zero angular separation which should equal the variance 
of the map. Fig. ^ shows the 2-point function obtained from a noise map that 
was obtained from the MAP algorithm. To a very good approximation, the 
noise in the sky map is uncorrelated from pixel to pixel. 



4 Power Spectrum Estimation from Full Sky Maps 

The most fundamental statistic to measure from the observed sky map is 
its angular power spectrum, C^, which measures the variance of the fluc- 
tuations over a range of angular scales. If the temperature fluctuations are 
gaussian, and the a priori probability of a given set of cosmological param- 
eters is uniform, then the power spectrum may be estimated by maximizing 
the multi-variate gaussian likelihood function 

, ^ exp(— im^ • C^^ • m) , ^ 

where m is a data vector (see below) and C is the covariance matrix of the 
data which has contributions from both the signal and the instrument noise, 
C = S + N. We can work in whatever basis is most convenient, in the pixel 
basis the data are the sky map pixel temperatures while in the spherical 
harmonic basis the data are the ai„i coefficients of the map. In the former 
basis the noise covariance is nearly diagonal, while in the latter, the signal 
covariance is 

Pixel basis: Spherical harmonic basis: 

Ti m ^ aim 

S ^ E£ ^^CePeicose,,) S ^ diag(C2, C2, . . . C3, C3, . . .) 
N afSij N iV(fm)(^m)'(see below). 

In the case of MAP, the length of the data vector, N^ata, is of order 1 mil- 
lion, so it is necessary to find methods for evaluating C that do not require a 
full inversion of the covariance matrix C, which requires 0{N^^^^) operations. 
As with the map-making, our approach is fundamentally an iterative one that 
exploits the ability to find an approximate inverse C~^. The most important 
feature in the data that makes this possible is the fact that we cover the 
full sky and that the galaxy cut we impose on the data is predominantly 
azimuthally symmetric in galactic coordinates. Of secondary importance for 
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this pre-conditioner is that fact that our noise per pixel is not too strongly 
varying across the sky. We discuss the pre-conditioner in more detail below. 
We maximize the likelihood by solving 

Ji- = = • C-i • • C-i • m + tr(C-i • P^ / = -2 \nC 

■pi — dC_ ^ ' 

— dCe 

using a Newton- Raphson approach to root finding that generates an iterative 
estimate of the angular power spectrum at each step 

Ce -C, --l^Fu,— (12) 

I' 

where Fu' is the Fisher matrix 



In order to implement the solution in (^13) we need to be able to evaluate the 

df 

ace 



following components of Fw £X- quickly 



tr(C-i-P^) (14) 
tr(C-i-P^-C-i-P^'). 

We use the spherical harmonic basis in which the data vector consists of 
the aim coefficients of the map obtained by least squares fitting on the cut 
sky. The signal covariance is diagonal in this basis, while the noise matrix is 
obtained from the normal equations for the aim fit 



^ ^ilL){e,ny"■ie^^y - V{trn) (15) 

{Im)' 

where 



^^(l-m)(lniy = l^i = l^i ■ 

The sums are over the uncut pixels in the sky map, and we have used the 
fact that the noise is uncorrelated from pixel to pixel. 



4.1 Evaluation of C • m 

The term • m appears repeatedly in the evaluation of (p^. We compute 
this by solving C • z = (S -|- N) • z = m for z. A more numerically tractable 
system is obtained by multiplying both sides by S2 • N^^ 

(I • • S5) . S5 . z = • • m • y (17) 
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where y is the spherical harmonic transform of the map, defined in ([l^). Note 
that y can be quickly computed in any pixelization scheme, such as HEALPix, 
that has the property of having pixel centers that lie on rings of constant 
latitude with fixed longitude spacing so that fast FFT methods may be used 
in the transform ||ll[, ||l^. We then solve ( p7| ) using an iterative conjugate 
gradient method with a pre-conditioner for the matrix A = (I+S2 - N^^ - 82 ). 
We find the following block-diagonal form of A to be a good starting point 

~ fl + S^-N-^-Si \ , . 

\ diag(I + S^-N-i-S^)y' ^ ' 

where is an approximate block-diagonal form of the noise matrix dis- 
cussed below. The lower-right block of A occupies the high £ portion of the 
matrix where the signal to noise ratio S2 • N^^ • S2 is low, so a diagonal 
approximation is adequate. The upper-left block occupies the low £ portion 
of the matrix where the signal dominates the noise, so we need a better es- 
timate of N~^. In practice we find this split works well at £ — 512 for the 
estimated MAP noise levels. As for the approximate form of N~^, defined 
in ([l6|), note that the dominant off-diagonal contributions arise from the az- 
imuthally symmetric galaxy cut, which couples different £ modes, but not m 
modes. Thus is approximately block diagonal, with perturbations in- 
duced by the non-uniform sky coverage of MAP. We therefore use a block 
diagonal form of N^^ as the pre-conditioner 

^(£m)(fm)' = ^(^m)(£m)'^™™'- ^^^^ 

Using the pre-conditioner dl^) we find that the conjugate gradient solution 



of (17) converges in approximately six iterations and requires only minutes 



of processing as a single processor job on an SGI Origins 2000. 



4.2 Evaluation of tr{C-^ • P^) and F«/ 

There are two approaches to evaluating tr(C"^ • P^). The first is to employ 
the approximate form tr(C~"'^ • P^) using the preconditioner dlq). The second 
is to note that since (m • m^) = C, it follows that tr(C~^ • P^= (m-^ • • 
P^-C~^ -m). Thus we can use Monte Carlo maps with the requisite signal and 
noise contributions to evaluate tr(C~^ -P^). This is the most computationally 
intensive part of the power spectrum estimation, requiring 0{NjncNiterN^) 
operations, where Nmc is the number of Monte Carlo realizations used, and 
Niter is the number of iterations used in the conjugate gradient solution of 
(p7|). In practice, we find it most efficient to use the approximate form for 
the first few iterations of ( p^ ) then switch to the Monte Carlo form for the 
final few. 

The same considerations can be applied to the evaluation of the Fisher 
matrix, however it is important to note that the Newton-Raphson method 
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does not require an accurate 2nd derivative in order to converge to the correct 
solution. Our approach to estimating errors in the final spectra are discussed 
more fully in JT^ . The result of applying this power spectrum estimation to 
a simulated 90 GHz sky map is shown in Fig. ^j. 
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Fig. 5. Sample Ci spectrum recovered from a simulated 90 GHz sky map with 
MAP-like noise properties. This spectrum can be shown to be an unbiased, nearly 
minimum variance estimate of the true power spectrum [Q. 
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